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Several static boundary value problems involving anisotropic media have been solved 
and the results are presented here. The solutions have been applied to simple electrode 
configurations to obtain formulas for in situ determinations of the resistivities of anisotropic 
conducting media. The boundaries treated include infinite, semi-infinite, slab and spherical 
shapes. The media are anisotropic with a common resistivity in two orthogonal directions 
and a lower resistivity in the third orthogonal direction. 



1. Introduction 

In the follownig', field solutions and measurement 
equations suitable for finding the electrical resistivi- 
ties of anisotropic conducting media and analogous 
parameters are derived for several boundary con- 
figurations. The results have been applied to tlie 
measurement of the resistivities of living muscle 
tissue in situ but they may also be used for the 
measurement of the electrical properties of other 
anisotropic media, for example; the conductivities 
of earth formations and semiconductors and the 
dielectric constants of uniaxial crystals. The meas- 
ured values of muscle resistivity and descriptions of 
the physical apparatus will be reported elsewhere. 

The anisotropy of muscle is correlated with the 
directions of muscle fibers. These appear to consist , 
electrically speaking, of high resistivity tubular 
membranes surrouniknl within and without with low 
resistivity fluids. The membranes have a common 
orientation in a given muscde, yielding one direction 
of low resistivity parallel to the fibers and two 
directions of high resistivity transverse to the fibers. 
To maintain the muscle in its natural state during 
the measurements, the measuring electrodes were 
placed on the muscle surface rather than in the 
tissue interior. With one exception, the equations 
developed here reflect these conditions. 

When applied to the measurements of skeletal 
muscle in living animals the precision of the methods 
to be described were limited by the inhomogeneities 
of the tissue and the difficulties of obtaining exact 
alinement of the electrode structure with the axes of 
anisotropy. 

2. Potential Field of a Spherical Current 
Source in an Infinite Anisotropic Medium 

The potential field of a splierical current source, 
and its specialization to a point source, in a linear, 
homogeneous, infinite, anisotropic, conducting medi- 
um with sink at infinity, forms the basis for all the 
results to be derived. For anisotropic media having 

1 Present address, Syracuse University, Syracuse, N.Y. 



one direction of high conductivity and two directions 
of low conductivity, the solution to the stated prob- 
lem can be obtained by a procedure given by Smythe. ^ 
The electrostatic equation of continuity can be 
written in rectangular coordinates as 

-V.J==(r,^ + cr, -^ + (7.-^=0 (1) 

in which J is the curnMit (hMisity, y the potential and 
ax, (Ty, and az are the conductivities in the x, y, and z 
directions respectively. The origin of coordinates 
is taken at the center of the sphere, figure 1. In tlie 
medium of interest, <Tx=(yz — (yi and (Ty=Gh. a~ is 
defined as (Til(Xh< 1. Equation (1) becomes Laplace's 
equation in a new coorchnate system in which 
lJL=^ay as 

2 Snivthe, W. R., Static and Dynamic p:icctricity, 2d ed. (McGraw-Hill Book 
Co., 1950), p. 248. 




Figure 1. Coordinate relations for calculatingjthe potential 
of a spherical electrode in an anisotropic medium. 
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(2) 



The spherical electrode surface defined by x^-^y^ 
-\-2^=Rq^ becomes, in the x, ^, z coordinates, an 
oblate spheriod defined by 



R% a^Rl Rl 



= 1. 



(3) 



Smythe shows that the solution to eq (2) may be 
expressed as 



'■-^"j.° 



dv 






in which v is the confocal ellipsoid para,meter of the 
equation 



{a'+v) ' {¥+v) ' {a'+v) 



= 1. 



(5) 



At the electrode surface defined hy v=Q, eq (5) 
reduces to eq (3). Hence a=Ro and b = aRo. By 
changing the variables of integration in eq (4) to 
v^ = a^-\-v, the integrals can be evaluated by Dwight ^ 



192-11, giving 



V=Vo 



tan 



.1 ib'-a'+v') 



l/2-|«' 



tan 



1 ib'-a^+y'V '-i'' 



2-1 a 

Ja 



(6) 



which reduces, upon insertion of the limits, to 



V= 



Vo 



■sin-'(l-a2)i/2 



tan 



_i {a' -by- 
{b^+vyi^ 



(7) 



Solving eq (5) for v and substituting in eq (7), 
yields the potential outside the spherical electrode 
in the original coordinates; 



V-- 



Fo 



sin-i(l-a2) 



,2\l/2 



Xsin-1 



(2^)^ 



[p'^+aY+AYliaY^P'^-Ay+^aYAY'Y 



12 



(8) 



in which A=R^^ {l-a') and p'={x''+zY^'> Equa- 
tion (8) also applies to a hemispherical electrode in a 
semi-infinite conducting medium; i.e., s>0 an 
insulator in figure (1). 

The point source solution is of particular interest 
and is found from eq (8) by assuming y and p' very 
large as 



V=^ 



sm-\l-a'y'' {p'^+aY) 



1/2 



Equation (9) can also be written in terms of the total 
current, /, leaving the electrode by integrating the 
normal component of J over a closed surface. By 

3 Dwight, H. B., Tables of Integrals and Other Mathematical Data, (Mac- 
Millan Co., New York), 3d edition. 



expressing J in terms of V and integrating over a 
cube of edge 2w, symmetrically disposed about the 
origin, we have 



/= 



v.iAy 



Ay r r' f^ 

l_,2)l/2L-.^J_J_(^2 

dy dz 



dx dy 



sin ^( 

Jw n 

dxdz 



'+ur'+aYy 



(V 



-y^'+aV) 



,2.,,2\3/2 



]• 



The bracket reduces to ^iriaiany^^. For a semi- 
infinite medium, ai and (Th = where 0>O and the 
bracket becomes 2T((Tia,y''^. Thus, in the latter 
case, eq (9) becomes 



y=/[47rV,a,(p'2+ay)]-^/2. 



(10) 



The measurement of a^, ay, and o-^ of a semi- 
infinite anisotropic medium can conveniently be 
made by applying eq (10) to a simple electrode 
arrangement. Assume the semi-infinite insulating 
boundary (usually air), figure 2, is at the z = plane; 
i.e., the conductivities are zero for z>0. Medium 
(1) is assumed anisotropic such that crx=(Tz=cri and 
Oy=(Th. The electrodes are shown schematically in 
figure 2 and consist of four, equally spaced, approxi- 
mately spherical, conductors of small size. The 
dimensions of the electrodes must be such as to 
make eq (10) a good approximation to eq (8) at a 
distance, a, from the origin; and to eliminate for 
practical purposes, perturbations by the remaining 
conductors of the potential field of the source or 
sink. The potential diflference between points g 
and h due to electrode e can be found from eq (10). 
Noting that electrode / doubles this potential differ- 
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Figure 2. Arrangement of four point-electrodes for measuriiig 
the resistivities of an anisotropic, semi-infinite, conducting 
medium. 

2>0 is an insulator. 



218 



ence and exprossiiif*; the result in polar coordinates, 
see figure 2, the potential diU'orcMice, Vgn, ean be 



written 

y,,(</))=/|47rVzc7, (eos2 0H-a^' siii"^ 0) j-^^ 

Two measurements of Vgn, at (/>=() and 
give (Ti and 0-7^ from eq (11) as follows: 

ai = I[2TaVgn{7rl2)]-' 
cTn=<ri[V,n{T/2)IVgn(0)Y^ 



(H) 

-7r/2, 

(12) 
(13) 



3. Measurement of the Resistivities of an 
Anisotropic Conducting Slab 

It may happen that the anisotropic material of 
interest, 3<^0, figure 2, has a limited extension in 
the negative 2 direction; the shape being referred to 
as a 'slab.' When slab thickness, h, is less than twice 
the electrode spacing, a, the effects ^of the boundary 
ma}^ make eqs (12) and (13) poor approximations 
to the slab resistivities. 

If the boimdary at 2=—h is either a perfect con- 
ductor or an insulator, however, exact mathematical 
expressions for the resistivity of the slab as a func- 
tion of quantities measurable with the electrodes of 
figure 2 can be derived. To obtain these expressions, 
we find first the following potential fields: 

The potential on the surface of an anisotropic 
conducting slab arising from a point current source 
on the surface and sink at infinity. Tlie slab lias 
a high conductivity in the y direction and low 
conductivity in the x and 2 directions, is bounded 
at 2=0 by an insulator and at 2=—h by (a) an 
insulator and (b) a perfect conductor. 

Generalizations of the method of solution to include 
arbitrary field point and source locations in the 
slab will be ajDparent. 

Figure 3a represents a section of the given problem 
in the y2 plane; figure 3b, a system of unages located 
in an infinite homogeneous medium having the slab 
conductivities. The image system satisfies the 
boundary conditions at s==0 and z=—h and can 
therefore be used to obtain the field within the slab. 
As shown in figure 3b the images extend to infinity 
in the plus and minus 2 directions, spaced a distance 
2h apart. They are equal in strength to the source 
plus its image in the 2 = Q plane; i.e., 2/ and are of 
the same sign \i2=—h is an insulator and alternate 
in sign if 2 = —h is a perfect conductor. The 
potential on the slab surface can be written by 
inspection using the principle of superposition and 
eq (10) suitably modified for a source not at the 
origin. The resulting potential on the surface 
z=0, is 



+2Z:(± l)V+ay+{2nhyT'''']- (14) 
«=1 J 



r = or c>o 



-^ 



^z--h 



(a) 




FicuRE 3. (a) yz-plane section of an anisotropic slab of thick- 
ness, h, with current source at the origin. 

(b) System of images of strength, rb 21, which 
match boundary conditions at slab surface. 

Images are infinite in number and spaced a distance, 2h, apart. The minus 
signs appear at alternate images when z<— /j is a perfect conductor. 



'I'lie plus sign applies when the 2=—h boundary is 
an insulator and the minus sign when it is a perfect 
conductor, eacli term in the summation being an 
image contribution.^ 

Again considering the electrodes of figure 2, the 
potential dift'erences measured on the slab surface at 
0=0 {y=0) and 0=7r/2 {x = 0) can be expressed as 

-[l + inh/ayr'^'ij^ (15a) 
T r 

Vghix 

-[l + {nhlaa)T"')'\ (15b) 



=o=2^[l+2S(±l)"(2[l + (2nAM)^]-'/^ 



The conductivities ci and (r„ can conveniently be 
found from the measurements of Vg^- Tlius, taking 
the ratio of the equations in (15) we have 



V, 



gh:y = 0_ 



aB 



F, 



gh.x = 



B{a) 



(16) 



in which B and B{a) are the terms in square brackets 
in eqs (15a) and (15b) respectively. Plots of B{a) 
as a function of hjaa (note that B=B(a) if a=l) 

* Although eq (14) diverges when the plus sign is used, eqs (15) do not. 
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can be found in the geophysical literature.^ The 
graphs are reproduced here in figure 4 with minor 
modifications. Using eq (16) and figure 4 it is 
possible to obtain a quickly by successive approxi- 
mations. As is apparent from figure 4, when 
a/A>1.5, the asymptotic formula, 5=1.39 {ajh), can 
be used; and when aa<.25, B{a) can be approxi- 
mated as unity. Once a={(Jil(jn)^'^ is found, eqs 
(15), together with figure 4 can be used to find 
a I and a^. 

4. Measurement of the Resistivities of an 
Anisotropic Conducting Hemisphere 

Should the anisotropic material have a spherical 
boundary, it is possible to obtain measurement 
equations in some instances. In particular, closed 
form expressions can be derived for the potential 
distribution in the herdisphere of figure 5. The 
source is a point current electrode arbitrarih^ located 
on the sphere surface and the sink a perfectly con- 
ducting meridian plane. The hemisphere has a 
high conductivity, an, in the <^ direction and low 
conductivity, cr^ in the R and 6 directions. Approxi- 
mate expressions, based on a four-electrode measure- 
ment, can be derived for an and o-/ from the potential 
field. 

This problem arose in measuring the resistivities 
of heart muscle. The muscle fibers were assumed 
to be wrapped around an axis in a manner roughly 
analogous to string wound into a ball. It should be 
noted that the field expression is unchanged if the 
hemisphere is replaced by a sphere and the con- 
ducting plane electrode b}" a point sink at the location 
of an image of the source in the plane. 

As shown in figure 5, the point electrode is on the 
surface at R=Cj d=di, <j)=<i)i and the plane electrode 
is at 4>=±Tr/2. The equation of continuity can, 
with an appropriate coordinate transformation, be 
expressed as Laplace's equation in spherical coordi- 
nates. The result is given in eq (17) in which 
<l)^ = a(j). 

V72T/ 1 ^ ^02 ^V\ . 1 ^ / • ^^V\ 



1 



dW 



--0. (17) 



' R' sin2 d 50'^ 

The boundary value problem in the new (jB, d, <i>') 
coordinate system is that of a point current source 
.at R=C, 6=dij <i>' — a(i)i on a wedge-shaped section 
of a sphere, and a sink at constant-potential inter- 
secting (internal angle ax) planes, figure 6. 

An image solution for the potential field in the 
wedge exists if 1/a is an integer; i.e., if the internal 
wedge angle is a submultiple of tt. Sources and 
sinks placed on the surface of a sphere as shown in 
figure 7 will make the planar wedge surfaces equi- 
potentials and satisfy the remaining boundary con- 
ditions as well. More general results arising from 



5 Hummel, J. N., AIME Tech. Publ. No. 418 and Heiland, C. A., Geophysical 
Exploration, (Prentice-Hall, New York, 1940), p. 716. 
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Figure 4. Graphs used to solve fc7- a in eq 16 (after Hummel). 




Figure 5. Anisotropic conducting hemisphere on a perfectly 
conducting meridian plane. 

Electrodes are on plane and at pomt, P, on the sphere surface. The hemisphere 
has a high conductivity, ah, in the 4> direction and equal but lower conductivities, 
<Ti, in the d and F directions. C is the radius of the hemisphere. 




Figure 6. 



Hemisphere of figure 5 transformed to new coordi- 
nates in which 0' = a<^. 
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a source inside^ the wcnlgo cnu bo ohtniiKMl hy similar 
reason iriti'. 

The poleiilial due to a source-sink pjiii" nihil rarily 
located in a sphere has been given by l^'iniik.'' His 
result, specialized to the case in which the source 
and sink as well as tlie field point are on the s])here 
surface and extended by superposition to q soui'ce- 
siidv pairs, is the desired etjuation for the potential 
arising from the system of images shown in figure 7. 
Til is extended result, which gives the potential on 
the wedge, is presented as eq (18). 

V=kp^±Q^-~ln [r'.+Cd-cos^k'd])- (18) 

O^the radius of the sphere. *^- 

fi^the chord distance from the ith source (or sink) 
point to the field point in the (li, d, </>') 
coordinates. ^a^^^"^ 

^-C[2(l-cos^;]^/2^ 
^i=the angle formed at the origin by the ra(hus 
vector to the field point (C, 6, </>') and the 
ith source (or sinlv) point {C, 6i, <i>'i). 
=cos~^ [cos 6 cos ^i+sin 6 sin di cos (</>' — </)')]. 
/:=l/27r(o-fO-;j)^^^, as can be shown by comparing 
eq (10) with eq (18) as C >oo. 

The potential field in the original coordinates is 
found bj substituting a<j) for </>' and a(j>i for </>■ in 
eq (18). 

When the conditions described below are satisfied, 
a convenient formula for determining approximate 
values of the conductivities ai and an, applies to the 
electrode configuration of figure 8. In this case, elec- 
trode e and electrode / each give rise to a set of 
images. The electrodes are spaced along the arc of 
a circle of radius, C, and the potential-measuring 
electrodes are at distances a and 2a from the current 



5 Frank, E., J. of Appl. Phys. 23, 1225, 1952. 



'"^1 + 




FiGXRE 7. Latitudinal plane section of figure 6 shoicing image 
construction for q = 3. 

The current, /, is shown entering at the point, P. 'I'lic radius of the circular 
section is, C sin^i. 



electrodes. If the source is located far from the 
constant-potential sink plane, figure 5; i.e., in a 
region defined by — 7r/4<(/)<7r/4, 37r/4<^<7r/4, and 
the electrode spacing a is small compared to the 
sphere radius, the contribution of the images to the 
potential difference measured by the electrodes of 
figure 8 can be reduced to tlie order of one percent. 
From (18) the potential differences, neglecting im- 
ages, measured by electrodes g and h, figure 8, due 
to the current electrodes e and /, reduce to conven- 
ient expressions at two part icular (^hu^t roch^ alinements ; 
i.e., when the electrodes He ah)iig a line of constant 4> 
or constant 6. In fact 



V. 



./.constant </> ^^(^^^^1/2 

x[(2[l-cos(^,-^,)l)-^/^^-(2fl-cos(^,-0])-^/' 



1 [2 -2 cos (^,-^,)]^/^+[l-cos {Be -e^) 
■^2 ^'' |2-2 COS (^,-^,)J^^^+[l-cos {d. 



)L 



(19) 



Here we have used de for the polar angle* of the 
source point, electroch* e] and 6^ arid On foi* the polar 
angles of the field ])oints, electrock^s g and h respec- 
tively. It is clear from the form of eq (19) and the 
angular relations sliown in figure 8 that the factor 
of two used in eq (19) correctly accounts for the 
effects of the shik electrode,/. Noting that the de- 
nominators of the first two tei'ins in bra.ckets on the 
right are chord lengtlis, a and 2a respectively, eq (19) 
becomes 



If/?: constant 



^ I (1 + 1 i„ 2(l+a/g) \ 



(20) 

In a similar wa}', wfien the elect roch's are aUned along 
a latitude line, and if the sines of angles a^ and 0:2, 
figure 8, subtended at the sphere center by the cur- 
rent and potential electrodes are approximated by 
the angles, w^e liave 



V, 



g;?: constant 6 



(-. 



1 



1 



. 2(l+aa/(?) \ 
"7r((r,<Tj"- \2aa/2C (l+ao,/2C)/ 



(21) 




Figure 8. Electrode arrangement for measuring the resistivi- 
ties of the hemisphere of figure 5 using eqs {20) and {21). 
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In the ratio of eq (20) and eq (21), a is the only 
unknown and may be found by successive approxi- 
mations; after which ai and an can be obtained from 
either eq (20) or eq (21). 

5. Resistance Between Two Small Spheres in 
an Anisotropic Medium 

An ^ ^apparent resistivity" of an anisotropic me- 
dium has been obtained by comparing the resistance 
measured between two small spherical electrodes in 
an anisotropic medium with the resistance between 
similar electrodes in an isotropic medium of known 
resistivity.^ When the two resistances are equal, the 
apparent resistivity is defined to be the known iso- 
tropic resistivity. The apparent resistivity, in addi- 
tion to being a function of ai and a^^ is a function of 
the orientation of the axis connecting the spheres. 
An approximate formula for this functional depend- 
ence will be derived. The coordinate relations are 
shown in figure 9. A first approxunation to the re- 
sistance, 1^, between two spheres is well known to be^ 



Q=2[V{B,)-V{&)]II 



(22) 



in which V{Rv) is the potential on the surface of the 
sphere of radius, Eq, centered at the origin and V{A) 
is the potential due to the sphere at the origin at the 
center of the second sphere (at radius vector d), 
neglecting the perturbing effect of the second sphere. 
Rewriting eq (8) in polar coordinates with x^-\-y^-{-z^ 
=(P and the constant Fo/sin~^(l — a-)^^^ as //47r 
{Gi(TnA)^''^j we have 



F(d) = 



/ 



.ww,..-'"-^^^^^r^' <^') 



^ Rchwan, U. P., and Kav, C. F., Circulation Research 4, 664 (1956). 
" Smytlic, W. R., op. cit., p. 238. 




Figure 9. Coordinate relations for evaluating approximate 
resistance between two spheres of radii, Ro, and center-to- 
center distance, d. 



in which 



-|-('ri-(l-a2) sm-^(^shr^^-(^^^Y (l-a^)]' 
+^a'^W4>^hv6i^) (l-a^)j 



y(iZo) is a special case of eq (23) obtained by letting 
d=Ro with the result 



V{Ro)= 



/sin-i (l-a^: 



^T{ffna,Ay 



(24) 



Comparing eq (22) for an isotropic medium in which 
Q> = (2Ta)~\llEo—l/d), with the equivalent isotropic 
expression obtained from eqs (22), (23), and (24), 
the apparent resistivity, pa, is given by 



sin ^ (1— o:^)^^^— sin 



Kd 



Pa= 



lUWn<Ji{l~a')f'\\IR,-\ld] 



(25) 



The two extreme values of K occur when the elec- 
trode axis is perpendicular and parallel to the y axis. 
In these cases, we have 



K 



= 0- 



= (2) 



1/2 



K,..„2A-.n=\2a^+2 (jfj (!-«')]'''■ 



(26) 



Insertion of eq (26) in eq (25) gives the extreme 
values of p^, hence some estimate of the mean of pa 
over all angles. The latter quantity can be ob- 
tained more precisely by nmuerical integration of eq 
(25) over d and </>. When Rold<C<^a, the terms in 
K containing Ro/d can be set to zero and the calcu- 
lation of the mean considerably simplified. Finally, 
by taking the ratio of eq (22) at the extreme values, 
an expression containing a as the only unknown can 
be obtained and from this ai and a^ determined. 
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